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Abstract: In this report we propose a MultiObjective (MO) performance 
evaluation framework for wireless ad hoc networks where criteria such as ca- 
pacity, robustness, energy and delay are optimized concurrently. Within such 
a framework, we can determine both the Pareto-optimal performance bounds 
and the networking parameters that provide these bounds. The originality of 
this approach is that it accounts for the inherent broadcast properties of the 
transmission and finely models the interference distribution. In the proposed 
model, the network performance can be optimized when several flows (source- 
destination transmissions) exist. One benefit of our approach is that the com- 
plexity does not grow with the number of flows. The other major contribution 
of this paper is the new analytical formulation of the performance metrics. It 
relies on a matrix representation of the constraints imposed by the interference- 
limited and broadcast wireless channel. Because of the similarity of this matrix 
with a Markovian transition matrix, we can exploit classical results from Markov 
chains theory to derive steady state performance metrics relative to capacity, 
robustness, energy and delay. Another very interesting feature of these new 
metrics is that the Pareto-optimal solutions related to them provide a tight 
bound on capacity, robustness, energy and delay. 
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De la prise en compte des interferences et de la 
nature broadcast du canal radio pour revaluation 
de performance des reseaux sans-fils 

Resume : Ce rapport de recherche presente un modele multiobjectif devaluation 
des performances d'un reseau ad hoc sans-fils. Les criteres tels que la capac- 
ity, la robustesse, l'energie et le delais de transmission sont optimises conjointe- 
ment. Grace a ce modele il est possible de determiner une borne en performance 
Pareto-optimale et les differents parametres du reseau qui permettent d'obtenir 
de telles bornes. L'originalite de cette approche reside dans le fait qu'elle modele 
finement la nature broadcast intrinseque du canal radio et prend en compte de 
fa con precise la distribution des interferences dans le reseau. II est ici possible 
d'optimiser les decisions de routage et d'ordonnancement des paquets dans un 
reseau quand plusieurs flots sont transmis dans le reseau. La complexity de cal- 
cul de l'ensemble des solutions Pareto-optimale n'augmente |pas avec le nombre 
de flots presents dans le reseau. La contribution majeure de ces travaux est la 
presentation d'un nouvelle formulation analytique des mesures de performance. 
Cette formulation se base sur une representation matricielle des contraintes im- 
posees par la nature broadcast du canal et la distribution des interferences dans 
le reseau. Du fait de la similarity de cette matrice avec une matrice Markovi- 
enne, il est possible d'exploiter des resultats classiques de la theorie des chaines 
de Markov pour deduire des metriques de capacite, robustesse, energie et delais, 
le tout pour un regime permanent. 

Mots-cles : reseaux ad hoc sans-fil, optimisation multi-objectifs, evaluation 
de performances 
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1 Introduction 

Wireless ad hoc and sensor networks are many times operating in difficult en- 
vironments and require several performance criteria to be satisfied, related to 
timely, reliable, and secure information transfer. Routing and resource alloca- 
tion protocols are key elements for ensuring the information transfer across the 
network. To better understand the capabilities of such protocols for a given 
network topology, it would be very helpful to know the bounds that can be 
achieved with respect to multiple performance criteria. These bounds illustrate 
the interdependence between multiple performance metrics and can capture the 
tradeoffs between the various operating points for the network. As a conse- 
quence, defining a unified multi-objective design framework capable of captur- 
ing the above mentioned tradeoffs for various possible operating points for the 
network becomes of premier importance. 

However, globally optimizing capacity for two-dimensional networks where 
routing and resource allocation (i.e. frequency, time or power assignment) are 
performed concurrently is a very hard problem that has triggered a comprehen- 
sive research effort under various conditions since the seminal work of Gupta 
and Kumar pQ. Several works provided ways for increasing the asymptotic ca- 
pacity bound of 0(n log(n)) [TJ[2] by for instance accounting for mobility [3]. 
In these works two important assumptions are made: unicast communications 
and threshold-based interference models. In a unicast communication, packets 
are sent to a specific receiver on a multi-hop chain. In a threshold-based model, 
nodes interfere with other nodes within a fixed range while beyond that range, 
no node is interfered. These results have been extended in [J] by accounting 
for realistic additive interference or in [S] by considering multi-user detection 
techniques. Toumpis et al. [5] have also proposed an interesting model to derive 
the capacity region by properly scheduling transmissions and hence accounting 
for a temporal multiplexing directly in a 2 dimensional network model. While 
the model is appealing, it is not scalable with respect to the network size and 
the number of flows transmitted. 

The broadcast nature of the wireless channel is mostly assumed as the main 
source of interference in these works and hence negatively impacts the network 
performance. However, this broadcast property can serve the capacity if the 
receiving nodes are able to cooperate and coordinate to transmit the data flow 
[7] . It has been shown to be beneficial in the context of opportunistic routing [5] 
as well. 

In this paper we consider the broadcast nature of the channel as a possible 
way to increase the capacity, and we propose realistic interference models. This 
work does not aim to assess the performance of a specific routing protocol, but 
more generally to extract the boundary of the feasible performance region of the 
network. Further, we consider that capacity is not the only interesting objective. 
More specifically, it is the trade-off between different performance criteria that 
should be considered to properly select a transmission strategy in the network. 
For instance, increasing the number of parallel paths improves capacity but 
at the expense on an increased energy due to the multiplication of packets 
traveling in the network. Thus, various criteria related to transmission delay [3], 
energy consumption [10] or fairness [11] should be considered in addition to the 
main design goal of reliable information transmission. As a consequence, the 
assessment of networking protocols usually relies on various criteria which may 

RR n° 7379 



4 



Jaffres-Runser, Gorce & Comaniciu 



be evaluated analytically or through network simulations. In this work, we 
concentrate on three criteria: capacity, energy and latency. As we will see in 
Section [SI robustness is related to capacity with a redundancy factor. 

We consider a network as a set of nodes comprised of relays, sources and 
receivers. In this work, we do not look for the optimal single path between a 
source and a destination, but for a set of relays deciding to forward packets or 
not. This approach is more closely related to a diffusion mechanisms than to 
a strict single path routing protocol. However, our more general formulation 
encompasses single path routing solutions. This concept has been inspired by 
the pioneering works on optic and electrostatic inspired routing |12H14) . How- 
ever, our work is not directly related to a physics inspired phenomenon. We 
provide here a more practical framework that relies on a time and spatial mul- 
tiplexed approach: we model the decision of each relay to transmit, receive or 
just sleep with a given probability associated to a time slot. These forwarding 
probabilities define a steady state of the network. Then, we propose an ana- 
lytical framework to obtain an estimate of overall capacity, latency and energy 
depending on the forwarding probabilities associated to each node and time slot. 
Since interference is modeled for a network steady state, the results obtained 
scale with respect to the number of flows being transmitted in the network. 

The paper is organized as follows. Section [2] gives some preliminary nota- 
tions. The network model is detailed in Section [3] and Section 0] describes the 
MO optimization problem. The steady-state performance evaluation tools are 
described in Section [3] and Section H2 concludes the paper. 

2 Preliminaries 

2.1 Notations 

Throughout this paper, upper case letters (e.g., X, Y, Z) usually denote vectors 
or matrices. Lower case letters (e.g., x, y, z) denote probability values. The 
transpose of a vector or a matrix is shown by F'. Subscripts denote the indeces 
of nodes of the network. For instance, Fi could denote a vector relative to node i 
and pij a probability related to both nodes i and j. Superscripts are introduced 
to specify the index of a time slot. For instance, x uv could denote a probability 
value relative to time slot u and v. Sets are denoted by calligraphic alphabets 
(e.g., A, B, C) and the cardinality of set X is denoted by \X\. The complement 
of a set A is shown by A c . 

A flow of symbols coming out of a node i on time slot u is denoted by /" 
and a flow of symbols coming into node i from node k in time slot v is denoted 
by fy. Table Q] summarizes our notations. 

2.2 Definitions of graph theory 

In this part, we briefly review the concepts and definitions from graph theory 
considered in this paper [15]. A complete graph Q = (V,£) has vertex set V and 
edge set £ C V x V. Without loss of generality, let V = {1, 2, ... , |V|}. 

We assume that the graph is finite, i.e., |V| < oo. For each node i g V, 
and M\ are the set of edges leaving from and the set of edges going into v, 
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Table 1: Some important notations in this paper 



v 

V 


uCU Ul IlULLCo 


p 


nf Ari cvac 
OCL (Jl CLINGS 




nf cniivpp nnriAC 


V 


Set of destination nodes 


K 


Set of possible relay nodes 


N 


Number of possible relay nodes 


T 


Set of time slots composing one frame 




Set of incoming edges at node i on time slot u 


$r? 


Set of outcoming edges at node i on time slot u 




Transmission rate for node i in time slot u 


A 


Discrete set of possible values of rf 


Pi 


Channel probability on edge (i, j) for slot u 


uv 


Forwarding probability for a packet sent by node i in 


time slot u to be forwarded by node j in time slot v 


M u 


Set of active transmission in time slot u 


M 


Set of all actives transmission 


fa, Id, /e 


Redundancy, delay and energy criteria, resp. 



respectively. Formally 

A* = {(i,j)\(i,j)€S} 

A complete graph /Ciyi of | V| vertices is a graph having the maximum number 
of edges. 

3 Network model 

We model the wireless ad hoc network by a finite complete graph /Ciyi for two 
reasons: modeling the broadcast nature of the wireless channel and finely ac- 
counting for interference. Transmission in the network is time multiplexed and 
a frame of \T\ time slots is repeated in each epoch s. Each edge (i,j) € £ for a 
particular time slot u represents an interference-limited channel which is mod- 
eled by the probability of a symbol or a packet to be correctly transmitted. This 
probability is referred to as the channel probability in the following. It models 
interference as an additive noise and is computed considering the distribution 
of the bit error rates (BER) or the packet error rates (PER) as shown here- 
after. As a consequence, there are \T\ orthogonal interference-limited channels 
for each edge € £ as illustrated on Fig.[T] Each channel is assumed to be 
in a hald-duplex mode, i.e. a node cannot transmit and receive a packet at the 
same time. 

A set of sources O and destinations T> is defined. A flow of packets is trans- 
mitted between any source-destination combination ({Oi, i € [1, .., |C|]}, {Dj,j € 
[1, .., |C|]}). We make the assumption that source and destination nodes do not 
relay the information. As a consequence, the network we are modeling is com- 
posed of a set of relay nodes 1Z = V — O — V. In the following, we consider that 
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the number of relays in the network is N ~ \1Z\. We also assume that a relay 
can not differentiate packets. As a consequence, all packets are treated as being 
unique by a relay. 

A transmission is defined as the couple (i,u) G V x T and represents the 
fact that node i is transmitting in a time slot it. 

3.1 Transmission rate 

We consider that a node j transmits a flow of symbols in time slot v. The 
symbols are the realization of a random variable X~j chosen from an alphabet 
X [X = {0,1} for instance). With this definition, we consider that a node 
j transmits the same symbol in a time slot v on all its outgoing edges Afj . 
This constraint incorporates the broadcast property of wireless communications. 
Time multiplexing is here accounted for in our model since a node can transmit 
different symbols on each available time slot. 

A flow of symbols coming into node j from node i on time slot it is modeled 
as a random variable Y£j. The symbols are transmitted by node i on the edge 
(i,j) G £ in time slot it and consequently, experience the probability pfj of 
being received successfully in j. Depending on the interference temperature of 
the network, symbols are received successfully or not at j. 

Having this, let Yj = {Yij,(i,j) G be the random variable giving all 
the symbols that can be received on the incoming channels of node j for all 



time slots where Yy = {Y$,(i,j) G A/>}. Let Xj = {XV, j G Aff] be the 
random variable giving all the symbols that can be transmitted on the outgoing 
channels of node j for all the time slots. The relation between the outgoing 



random variables Xj and the incoming random variables Yj defines a coding 
scheme for the network. 

A node has one main decision to take upon receiving a symbol or a packet: 
whether it should process it or not. If it decides to process it, the next steps are 
to 

• Decide on the coding for the symbol or packet to be sent. 

• Decide on which channel to transmit it. 

These decisions strongly influence the quality of the transmission and aim 
at mitigating the transmission errors due to fading and interference. Depending 
on these decisions, the transmission rate of a node i on a channel u varies. For 
instance, if a node decides to drop one packet out of two received on a same 
time slot it, the transmission rate on channel it would become half the rate at 
which it received packets on the same time slot. For instance, it is possible that 
a node transmits one packet every two received packets on time slot u because 
it is applying a coding scheme for which two packets are combined into a single 
transmitted packet. In both cases, the node is adapting the transmission rate 
on each channel to combat errors. 

Based on this statement, we characterize the behavior of a node i by the 
rate at which it is transmitting in each time slot it. The transmission rate of 
node i in time slot u is denoted by r". Having this, a vector of transmission 
rates for each time slot can be defined 




\T\ 



, r? G [0,1 
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Figure 1: Relay network model: transmission rate and channel probabilities 



The transmission rate can be interpreted as the percentage of time a node is 
transmitting packets in a given time slot. For instance, if a node i transmits at 
a rate of r" = 0.5 in time slot u, it transmits a packet every two time slots, if 
one time slot permits to transmit exactly one packet. We assume that when a 
node is not transmitting on a channel, which happens with probability 1 — t", 
it is listening to packets. This is to comply with the half duplex assumption. In 
other words, the listening rate is equal to 1 — rf . The model could be extended 
to also account for the time a node may be sleeping, but this option has been 
disregarded so far because it would increase the dimension of our search space. 
However, it would be worth investigating such an option when minimal energy 
consumption prevails as it is the case for wireless sensor networks. The proposed 
network model is consequently defined in a steady state mode, where we know 
which nodes are transmitting or not packets, continuously. 

The average rate at which all the symbols are coming into node j is given 

by 

^7 = E E «■ 

uer (i,])eJf} 

for pfj the probability to receive a symbol correctly. The rate at which symbols 
are being transmitted by node j is given by: 

^ = E-J 

Let 

be the matrix of the transmission rates assigned to all the nodes of the network 
for all time slots. A particular instance of r belongs to the set of possible 
transmission rate matrices shown by T. r is feasible if Properties 1 and 2 hold 
for each node: 

Property 1: Flow conservation. The rate of all outgoing flows is lower or 
equal to the rate of all incoming flows, i.e. 

7]<r 3 (i) 

Property 2: Half duplex. A node j is able to receive a message on a time 
slot u the proportion of time it is not transmitting on that same time slot. As 
a consequence, the proportion of packets being actually received on time slot u 
is given by 

^(1-t/)+t? < 1, VueT (2) 
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where rt? = Y], , tt^t^p^ stands for the incoming rate in time slot u. The 

conditions of ^ are always true for the definition of the transmission and 
listening rates. However, this constraint becomes more meaningful if a sleeping 
rate variable is defined for a node. 

Now that we have defined the transmission rate matrix, we can define the 
set of active transmissions M. corresponding to r. A transmission (i, u) 6 V x T 
is said to be active if r" > 0. As a consequence, the set M. is given by: 

M = {(*,«) G V x T | t? > 0} 

Similarly, the set M. u refers to the set of transmission being active in a same 
time slot u. We have M u = {(i,!i)eVx T|t" > 0, v = u} 



3.2 Channel probability 

Knowing r, it is straightforward to derive the interference temperature of the 
network. We model interference as an additive noise computed at the end of 
each link (i, j) S £ of the network. Let pfj be the probability for a symbol trans- 
mitted by node i to arrive successfully at node j in time slot u. We denote p^- 
as the channel probability. It is a function of the statistical distribution of the 
Signal to Noise and Interference Ratio (SINR) at the location of the destination 
node j. This quantity can either be defined using a Packet Error Rate (PER) 
or a Bit Error Rate (BER) depending if transmissions on the network are pack- 
etized or not. In the following we consider that transmissions are packetized. 
Before giving the exact expression of the channel probability, a few preliminary 
definitions and notations are given hereafter: 



Pathloss attenuation factor and transmission power reflects the at- 
tenuation due to propagation effects between nodes i and j. In our simulations, 
the simple isotropic propagation model is considered. We consider that all the 
nodes use the same transmission power denoted as Pp. 



Interference Since we consider time-multiplexed channels, interference only 
occurs between transmissions using the same time slot. Let Zy be the power 
of the interference on the link € £ on time slot u and computed at node 
j. If we denote by a set of nodes of the network interfering at node j (not 
including i) in time slot u, interference power is defined by: 

T ij = £ P t ■ «W(fe)i for Id(fc) ^ * (3) 
where Id{k) gives the number of the node interfering at j. 



SINR The SINR between any two nodes i and j in resource u is given by the 
following equation: 

where ijy is the interference power on the link and Nq the noise power density. 
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Packet error rate (PER) For a specific value of SINR 7, the packet error 
rate PER can be computed according to: 

PER(~f) = 1 - [1 - BER(j)] Nb (5) 

where Nb is the number of bits of a data packet and BER(j) is the bit error rate 
for the specified SINR per bit 7 which depends on the physical layer technology 
and the statistics of the channel. Results are given for an AWGN channel and a 
BPSK modulation without coding where BER(j) = Q U/2y) = 0.5 * erfe( N /7). 



Interfering sets A node i is said to be active in the network if rf > 0. We 
recall that M. u designates the set of active transmissions in time slot u. We give 
now a more precise definition of an interfering set. An interfering set X^(K) on 
a link (i, j) in time slot u is any subset of the set M. u — {i} made of K elements, 
K = {1, .., \M u -{i}\}. If M = \M u -{i}\, let refer to the set of all possible 
interfering sets of cardinality \Cfj \ — Y1k=i Cif) + 1- 

Equation © details the derivation of the channel probability pf^ as the 
average of the PER experienced for all possible interfering sets / £ on time 
slot u referred to as PERf. 

where PERi = PER(ji), I <S Cfj, where 7; is the SINR experienced on the link 
in time slot u if the nodes of interfering set I are active. 
P; is the probability for the interfering set I to be active and create interfer- 
ence on the link on time slot u. More specifically, it is the probability that 
the nodes of the interfering set I are transmitting concurrently and the others 
are not as specified in the following: 

K M-K-\ 

p < = n r ^ n o- - o (7) 

fe=l m=l 

In (|7J|, IlfcLi T k gi ves the probability that the K active nodes of the interfering 
set I are transmitting and Ilm^i*' -1 0- — T m) t ne probability that the M — K — 1 
other active nodes are not. 

Similarly to the vector of transmission rates, we define the vector of channel 
probabilities for a link € £'■ 



1 \T\ 
Pij ■■■ Pij 



,p£e[o,i] 



Let Pj be the matrix giving all incoming channel probabilities at node j 

t 



Pi 



pt ... pt 



and P the matrix of all channel probabilities in the network 

P=[ Pi ■■■ Pn ] 
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3.3 Forwarding and scheduling 

Now, lets introduce the forwarding and scheduling decisions of the nodes relative 
to the transmission of a packet. This decision is represented by the probability 
of a node j to transmit on time slot v a packet coming from node i on time 
slot u. We will refer in the rest of the text to the forwarding probability xfj '• 
For each node of the network, we can define a (./V.|T|)-by-|7~| matrix giving all 
the forwarding probabilities relative to any node j of the network as follows. 
This forwarding matrix Xj is given by: 

X 3 = [ X lj ' ' ' X Nj ] f 

where each matrix provides the scheduling probabilities of a flow of packets 
coming from node i on its output times slots, depending on the time slot the 
packets are received on. 







i\n\ I 
• x ij 


X^ — 










\ni\n\ 
■ x ij 



The matrix of forwarding probabilities is related to the matrix of transmis- 
sion rates r and the matrix of channel probabilities P with the following set of 
equations 

E E^-^-^K- = T h v o»e^ (8) 

where r"p^ is the probability that a packet sent by i on time slot u arrives and 
(1 — Tj) is the probability that node j is listening on channel v. These equations 
introduce strict constraints on the choices of the forwarding probabilities. 

The forwarding probabilities represent the decisions of the nodes to either 
(i) retransmit all the packets or symbols received or (ii) reduce the output rate 
by dropping or re-encoding them together. From now on, we will refer to the 
set of all forwarding probabilities of the complete network using a matrix X of 
size i\r.|T|-by-iV.|T| defined by: 

X = [Xx . . . X N ] , X e X (9) 
where X is the set of all possible matrix instances. 

4 MO optimization problem 

The scope of this section is now to take advantage of the previously described 
network model to derive a framework capable of extracting the set of Pareto- 
optimal transmission strategies with respect to various performance criteria (e.g. 
capacity, delay, energy. . . ) . A Pareto-optimal set is composed of all the non- 
dominated solutions of the MO problem with respect to the performance metrics 
considered. The definition of dominance is: 

Definition 1: A solution A dominates a solution B for a n— objective MO 
problem if A is at least as good as B for all the objectives and A is strictly better 
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fi 

Figure 2: Non-dominated solutions are highlighted with red circles for a 2- 
function minimization problem. 



than B for at least one objective. Mathematically, we have for a minimization 
problem: 

A>~ B = 

Vi G [l,n] : fi (A) < fi(B),3j G [l,n] : < 

Notation Ay- B means that A strictly dominates B. In the space of the evalua- 
tion functions, for the case of a minimization of all the optimization objectives, 
the set of Pareto-optimal solutions for an example 2-objective problem is illus- 
trated in Fig. [2] 



4.1 Solution definition 

A solution of our MO optimization problem is given by the set of all the forward- 
ing probabilities represented by matrix X G X defined in ^ . As a consequence, 
the set X becomes the problem search space. Let fc, fa, fo and /g be the per- 
formance criteria relative to capacity, robustness, delay and energy, respectively. 
Capacity and robustness are maximized, while energy and delay are minimized. 
The derivation of these metrics is provided in the next section. Our goal is to 
solve the following multiobjective optimization problem by finding the set of 
Pareto-optimal solutions X opt : 

X opt = {A G X | VA G Xopt, VS G X£ pt , A y B} (11) 

where X opt U X£ t = X . The dominance relation is defined in (|10p and adapted 
for criteria fcfii, Id and f E . 

The cardinality of the search space X as defined grows exponentially with N 
and \T\. However, considering the complete search space may be meaningless 
depending on the network connectivity. As a consequence, it is reasonable to 
define the search space for a maximum number of relays N max < N . That is 
why we define Xjss ma:c as the subset of X made of at most N max active relays 
and refer in that case to the A max -relay problem. In this case, we are looking 
for the Pareto-optimal combinations X of N max relays in the set 1Z. 

There is a whole set of constraints for a solution X G <-tjv max to be valid 
related to Properties 1 and 2. As a consequence, the solutions of Xj^ mai that 
do not follow these properties are dropped before being evaluated by the search 
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algorithm since they are considered as not feasible. We recall that Property 
1 ensures that the flow conservation constraint at each node holds and that 
Property 2 ensures that the half duplex constraint is met. 

4.2 Change of variable 

The derivation of the transition rate matrix r € T knowing a solution X € X 
is intractable. In a nutshell, to compute the transmission rates of a node j, 
X and the incoming transmission probabilities Pj have to be known yet the 
elements of Pj are a function of the transmission rates matrix, creating a circular 
dependency between the variables. 

Nevertheless, being able to model exactly the interference temperature of 
the network is appealing and to keep the same formalism, we propose a reverse 
approach where a solution of the optimization problem is defined by the set r 
of the transmission rates for all the nodes and time slots. From r, it is possible 
to derive the channel probabilities matrix P according to ([7]) since the activity 
of all nodes on each time slot is known. 

Only instances of r that meet the constraints relative to Property 1 and 2 
are further considered as valid. Now that we have a valid r, we can derive all the 
forwarding matrices X e X that verify the constraints of equation ([5]). There 
are constraints, each one constraining the choice of the x™ for all nodes 
and time slots of the network with respect to r. Let X T be the subset of X 
that verifies ([5]) with respect to the transmission rate matrix r. Each solution 
X € X T can be evaluated according to fc, fa, /d and Je- 

The MO optimization problem of (fTTj) stays unchanged, however, the way 
the search space is constructed has changed and permits to select in a first stage 
a subset of valid and interesting solutions X T . This feature is very interesting 
in terms of optimization. In fact, if it is possible to determine knowing t if 
the set of corresponding X T is promising or not, it is possible to apply pruning 
techniques that disregard non-promising t's instances. 

4.3 Problem complexity 

We recall that each transmission rate r" takes its values in the continuous closed 
set [0, 1]. However, to be able to use common MO optimization heuristics, we 
formulate the problem as a combinatorial MO optimization problem. Therefore, 
transmission rates take their values in a discretized set A = {Si, . . . , 5|a|}i $i G 
[0,1] of |A| values. 

The set T of possible r solutions has a cardinality of |T| = |/A| JV -I' ?- I . The 
set of feasible solutions is reduced by Property 1 but since the relation between 
r and P is complex, it is difficult to derive analytically the exact number of 
feasible solutions. We can just say that T is bounded by OdAI^ I 7 "'). 

The set X of possible forwarding probabilities matrices has a cardinality of 
\X\ = |A| 2Ar ( JV+ l C) l- 1 )'l r l. The subsets X T of feasible X variables is also reduced 
because of Property 1. Let Xf be the set of feasible X solutions. Xf is the union 
of the X T for all feasible r instances. Similarly to T, it is difficult to evaluate 
the size of Xf. Thus, X f = 0(\A\ n2 -^). 
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Table 2: Main notations in this section 



F~Z(s) 


Flow vector at time epoch s 


M 


Transition matrix 


Q 


Relaying matrix 


D 


Arrival matrix 


M F 


Fundamental matrix 


Sm 


Source matrix 



5 Steady state performance evaluation 

This section provides a framework for the definition of performance criteria for a 
particular solution X £ X of our MO optimization problem. It is an important 
contribution because it permits a fast performance derivation for end-to-end cri- 
teria such as capacity, robustness, delay and energy. Moreover, this framework 
also permits to reduce the complexity of the problem by pruning solutions where 
the transmission rate matrix r does not meet sufficient performance constraints. 
Our framework relies on the definition of a transition matrix which is composed 
of the probabilities for a flow of packets to be re-transmitted by the nodes of 
the network. Such a formulation is inspired by the theory of Markov chains. 
However, we want to stress that we do not model the network using a Markov 
chain. We define instead a transition matrix that has the properties needed to 
be able to re- use some results from the theory of Markov chains. Table [2] gives 
notations specific to this section. 



5.1 Flow vector 



Let F™(s) represent the probability of node i to transmit a packet in time slot 
u at time epoch s. We consider that at the beginning of a time epoch, a frame 
of |T| time slots starts and all nodes having a packet to send can transmit it 
in their respective time slots. A new epoch starts when all nodes have finished 
receiving their packets and are ready to send new ones if needed. Let F out (s) 
be a vectorial representation of these probabilities for each node and each time 
slot at a given time epoch s. F out (s) is referred to as the flow vector and is 
shown by: 

F^ t (s) = \ f\( s ) ... F^(s), Fb(s) 



where f\(s) 



^(-) 



The vector F D (s) stands for the flows received at the destinations nodes of 



set V. It is given by F D (s) = F D (s) 



F r 



where F D .(s) gives 



the flows received in each time slot by destination Di . 

The flow vector is a (N + |2?|).|T| matrix when all the possible transmissions 
are accounted for. However, if a relay i never transmits in a time slot u (i.e. 

t" =0), we have F"(s) =0. As a consequence, the size of the flow vector can 
be reduced to the size of the only transmissions that are active. Hence, the size 
of F out (s) can be reduced to a \M \ + \V\.\T\. 
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5.2 Transition matrix 

The propagation of the packets and the forwarding decisions of the nodes can 
be represented by a transition matrix M. Each element of the matrix M gives 
the probability for a packet transmitted by one relay in one of the time slots 
to be transmitted by any other relay in the network in one of the time slots. 
When it is applied to the flow vector of time epoch s. the new flow vector of 
time epoch (s + 1) is obtained. Formally we have: 



F out (s + 1) = F out (s) ■ M 
The transition matrix M has the following structure: 

M = 



\ Q 


D ' 





I 



For I = N.\T\ and m = \T)\\T\, Q is a non-zero l-by-l matrix, D is a non-zero 
l-by-m matrix. / is a m-by-m identity matrix and a m-by-Z zero matrix. Q is 
referred to as the relaying matrix and D as the arrival matrix. 

Transition matrix M has a similar canonical structure than a finite absorbing 
Markov chain [16] . Here, we have I transient states, i.e. relay nodes that are 
forwarding packets using active transmissions. The relaying matrix Q gives the 
probabilities for the packets sent at time epoch s to be retransmitted at time 
epoch s + 1 by the relays of the networks. We have m absorbing states where 
the destination nodes of D receive packets on each time slot and keep them. 
The identity matrix represents the fact that packets received by a destination 
are never forwarded, but absorbed. Matrix D is composed of the probabilities 
from going from one transient state to one absorbing state, i.e. the probability 
for packets to be received at the destinations when being transmitted from any 
relay or source in the network. 

The relaying matrix Q is structured as follows: 



Q = 



Q 12 
Q21 



'Nl 



Q 



N-1N 



*1N 
*2N 





is an |7~|-by-|T| zero matrix representing the fact that a node i does not 
forward a packet coming from itself. The matrix is a |7~|-by-|7~| matrix 
that gives the probabilities of node j to transmit a packet sent by node i for all 
possible combinations of time slot. Its structure is: 



(12) 



where Qf? is the probability for a node j to retransmit on channel v a packet 
that has been transmitted by node i on time slot u. With respect to our network 
model, it is equal to: 







■ e 1 


ij 




• Qir iiri . 



p«(i 
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It is important to note that the powers of matrix Q tend to if its elements 
verify the strict inequality < 1. For this property to hold, we have to ensure 
that pfj < 1. This is usually the case for a wireless transmission where perfect 
channels exist in very few unrealistic cases. 
The arrival matrix D is given by: 



D 



D iD . is a |T|-by-|T| diagonal matrix whose diagonal elements Df D give the 
probabilities for a packet transmitted by a node i in time slot u to arrive at 
destination Dj. Consequently, D" D is the probability of reception: 

DiDj = PiDj (1 — T L> 3 ) = PiDj 

since t d . — 0,Vw G T for any destination node Dj. 

The source matrix Sm is needed to compute the initial flow vector F out (l), 
naming the probabilities of all the nodes to transmit a packet sent by a source. 
Sm is a |0|-by-|T| matrix. An element Tg. gives the probability of the source i 
to transmit in time slot u. An initial flow matrix is obtained from Sm with 

F)Z(1) = S M M S 

where Ms is a \T\-by-(N + |2?|).|T| transition matrix defined as Ms = [Qs D$] 
with Qs and Ds the relaying and arrival matrix for the packets sent by the 
source, respectively. We have 

Qs = [ Q S i ■ ■ ■ Qsn ] 

and 



Ds = 



Dsd x ' ' ' D SD 



where Q S i follows the pattern given by (fT!2|) . Since destinations do not retrans- 
mit packets, Dsd { is a |7~|-by-|7~| diagonal matrix whose diagonal elements are 
D u SDi =P u SDi ,u€T.^ 

The flow matrix F src (l) created at epoch s = 1 is a \0\-by-(N + |£>|).|T| 
matrix of flows. The initial flow vector F out (l) relative to source i is given by 
the i th line of F src (l). A very interesting feature of the steady state performance 
evaluation and network model proposed herein is that the matrix Q and D can 
be applied to any type of source-destination communication (unicast, multicast, 
anycast, ...). Indeed, this is a direct consequence of the fact that the inherent 
broadcast and interference- limited properties of the radio channel are accounted 
for in the structure of Q. For instance, for a multicast transmission, the source 
matrix Sm is composed of a single source and the arrival matrix D is defined 
for \T>\ > 1 destinations. The matrix Q stays unchanged. Or when multiple 
sources are transmitting to a same destination (a sink for instance in a wireless 
sensor network), a source matrix of Sm is composed of several nodes. Each one 
of them is propagated in the network using the same Q to a unique destination 
modeled by D. 
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For clarity purposes, we have presented all matrices using all N possible 
relays. However, when considering a specific solution t, these matrices are 
derived only for the set M. of active transmissions. For instance, the size of the 
square matrix Q is reduced from I = N.\T\ to I = Ad. 

5.3 Fundamental matrix 

We really want to stress that the proposed transition matrix definition is not the 
definition of a Markov chain. Indeed, the sum of the probabilities on a line of 
M is higher than 1. However, having Q™ < 1 ensures that Q s tends to zero as 
the number of time epochs s tends to infinity. As a consequence, the following 
theorem still holds: 

Theorem 1: I — Q has an inverse, and 

oo 

(I-Q)- 1 = I + Q + Q 2 + --- = Y,Q S 

s=0 

PROOF The complete proof is omitted for conciseness purposes but can be 
found in [16] . Even though M is not a Markovian transition matrix, the only 
property needed for Theorem 1 to hold in our case is that Q k — J- as k — >• oo. 

q.e.d. 

The inverse of / — Q is defined as the fundamental matrix and denoted 
Mf = (I — Q)^ 1 ■ In our case, the fundamental matrix gives us the mean of the 
total number of transmissions that have been done until all packets be received 
at the destination. 

5.4 Optimization criteria 

The definition of our optimization criteria are directly derived from the fun- 
damental matrix. Indeed, if I — Q is invertible, Mp gives a measure of the 
performance of one solution X G X when the number of time epochs tends to 
infinity. The very interesting feature is that it also accounts for all the possi- 
ble cycles in the graph and consequently models the broadcast property of the 
wireless channel. The criteria here after are defined for one source-destination 
flow. They will be extended in future work to multiple flows which is possible 
since the complete framework can model concurrent flows on the network. 

Redundancy, capacity and robustness It is possible to derive for one 
flow the number of packets received at destination for one transmitted packet 
(\S M \ = 1) with 



f = F out (l)- 



' M F D ' 




' " 


/ 




1 



where [ 0^ 1^ ] is a vector of iV|T| zeros followed by |2?||7~| ones that accu- 
mulate all packets being received at the destinations into /. This model shows 
that a solution X £ X provides / packets at the destination for one sent packet. 
However, since relays are not able to discriminate packets, the set of / packets 
may be composed in the worst case of / copies of a packet by the source. What 
we can say here is that: 
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• if / > 1, there may be at most / different packets. Consequently, / 
provides both an upper bound on the network capacity and redundancy. Let 
fc be the bound on the real capacity of the network. It is given by 

f c = min(l,/) 

This bound can be reached if all the packets received at the source are orthogo- 
nal. This may be possible if an optimal network coding strategy is found where 
all transmitted packets are encoded combinations of the original information. 

• if / < 1, / can be interpreted as a robustness criterion giving the proba- 
bility to receive one message at destination. 

This criterion can be used to remove solutions that do not guaranty robust- 
ness. For instance, solutions having / < 1 can be disregarded if robustness 
wants to be guaranteed. 

Delay criterion Let p(H = h) be the probability for a transmission to be 
done in h hops. Thus 



p(H = h) = 



D 
I 



' Q 




' D ' 







I 



h = 2 
h>2 



We assume that a relay introduces a delay of 1 unit. Consequently, a h-hop 
transmission introduces a delay of h — 1 units. It is given by the infinite sum: 



Id+2Q 



-{h+l)Q 



h=l 
h , 







-1) 


■p(H = 


h) = 




' D ' 




' " 




I 




1 



Let V be equal to (Id + 2Q-\ \-{h + l)Q h + ••). We can show that V = {M F ) 2 

and consequently 



Id - F out {\) 



by factorizing Q: 



(M F ) 2 D - 




' " 




I 




1 




Mp in the definition of V 


3Q 2 + • • • + 


(h 


+ l)Q h + ■ 



We now have V = M F + Q ■ V. Consequently V = (M F f. 



q.e.d. 



Energy criterion The energy is measured by the number of packets being 
transmitted in the network. It can be simply derived by: 



fE - F out {\) ■ 



' M F ' 




' 1 " 











where 1 is a iV|T|-by-l vector of ones that sums all the packets sent by the 
relays that have participated in the transmission. 
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td PSD 

Figure 3: 1-Relay network 



5.5 Example: 1-relay network 

In this example, we consider a basic 1-relay network as the one depicted in Fig. [3] 
Interference is completely mitigated in this network if we use |T| = 2 times slots 
and S transmits in time slot 1 and relay R in time slot 2. As a consequence we 
have Sm = Wg 0] and tr = [0 r^l . The forwarding probabilities matrix is: 



.11 

l SR 



b SR 



,12 
SB 



^SR 



tI/{t s . Psr .{1-tD) 




where the values of the xfj are derived according to ([S])- The derivation of 
Kuti 1 ) = SmM s provides: 

KlO) = [ ^s-Psr-^ ~ r 2 R ).xf R t s . Psd ] 

Since we have only one relay, matrix Q does not exist and we have M = 
[Z)t /t] t . The criteria / equals: 



' D ' 




' 1 ' 


I 




1 



= F 0Ut {l) 



Prd 

1 
1 



where / is obtained for Tg = 1. The bound on capacity of this 1-relay network 
can be obtained for / = 1. Thus, we can derive the value of (r^)* that guaranties 
/ = 1 and verifies: 



PSD + (T 2 R Y 



■Prd 



1 



It shows that the optimal transmission rate of the relay with respect to capacity 
is the rate at which the transmissions on the link (R, D) compensate the losses 
on the link (S,D). However, this bound may not be achievable when the sum 
of the information coming into D and R is lower than 1. Typically, applying a 
cut-set between S and the set {R, D}, leads to a maximal flow / < 1— (1 ■ 
(1 — Psr)- Further, this capacity bound can only be achieved if the relay knows 
which packet have not been received by D directly from S or, more interestingly, 
by performing a network coding to introduce proper diversity in the flow. 



6 Conclusion and discussion 

The work presented in this paper has derived a flexible and powerful framework 
for evaluating the performance of a wireless ad hoc network with respect to 
several performance criteria. It has been designed to account for the broadcast 
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nature of wireless communications and for an accurate interference temperature 
characterization for the network. The proposed framework allows for the de- 
termination of Pareto-optimal solutions with respect to capacity, redundancy, 
delay and energy optimization criteria. The solutions defined herein represent a 
steady state of the network where all active relays are constantly transmitting 
data. 

The complexity of the proposed framework is not a function of the number 
or the structure of the source-destination flows and it is scalable with respect 
to the number of flows in the network. As a consequence, it can model various 
types of transmission: one-to-one, one-to-many, many-to-one or many-to-many. 

This flexibility comes at the price of an extended problem search space. How- 
ever, the structure of this search space is very interesting from the optimization 
point of view since there is a dependency between the set of feasible transmis- 
sion rate matrices T and the set of feasible forwarding probabilities matrices 
Xf. If it is possible to disregard a feasible value of r because it is not promising 
enough for the MO search, then the worst case complexity of the problem can 
be reduced from 0(\A\ n2 W) to 0(\A\ N W). 

The definition of the fundamental matrix permits to derive an upper bound 
on the capacity of the network. This upper bound is provided in terms of 
the average number of packets being received at the destination. With our 
framework, we can not differentiate packets at the destination and hence can 
not state whether a received packets is a copy of an already received packet 
of not. However, we state that using proper network coding techniques, it is 
possible to get an network effective capacity close enough from this bound. 

This paper has focused on presenting theoretically our MO performance 
evaluation framework. It is only illustrated by a basic example. In future work, 
we will have to derive proper MO optimization search algorithms |17| to effi- 
ciently compute the set of Pareto-optimal solutions. This task is challenging 
but as mentioned in the paper, it is often enough to address the A^„ a2; -relaying 
subproblem (with N max < N the maximum number of active relays in the net- 
work). Next, it will be of major interest to analyze the Pareto-optimal bounds 
and its corresponding solutions. The knowledge of such bounds will provide a 
powerful tool to assess the quality of distributed routing and resource allocation 
algorithms. Further, and this is of great interest in our opinion, we can take 
advantage of the steady state Pareto-optimal solutions to design distributed 
protocols that get close to our MO performance bound. 
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